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The statistics of energy levels of electrons in a random potential is considered in the critical energy 
window near the mobility edge. It is shown that the multifractality of critical wave functions results in 
the violation of the normalization sum rule in the thermodynamic limit and leads to the quasi- Poisson 
term ((5N) 2 } = a(N) in the level number variance. The sum rule deficiency a = n/2d is related to 
the multifractality exponent r\ = d — d*(2) in the d-dimensional space. 



I. INTRODUCTION 

The energy level statistics in complex quantum systems are the subject of an advance 
study since the pioneer works by Wigner, Dyson and Mehta that led to the development of 
the classical Random Matrix Theory (RMT) 1 \ The energy level distribution found in the 
framework of the RMT, known as the Wigner-Dyson (WD) statistics, posesses a remarkable 
property of universality: it depends only on the symmetry of the Hamiltonian with respect to 
the time-reversal transformation T. There are three symmetry classes of Gaussian ensembles 
of random matrices: orthogonal, unitary, and symplectic, labeled by /3 — 1, 2 or 4, respectively. 
Orthogonal and symplectic ensembles correspond to time-reversal invariant Hamiltonians with 
T 2 equals +1 or — 1, and the unitary ensemble corresponds to the case of broken time-reversal 
symmetry. These three random matrix ensembles turned out to describe spectral correlations in 
a great variety of complex quantum systems, from nuclei to disordered metals, whose classical 
counterparts exhibit a chaotic behavior. 

Other universal statistics which apply to describe the spectral correlations in systems with 
non-chaotic classical motion, are the well known Poisson statistics. The principal difference 
between WD and Poisson statistics is that there is a phenomenon of level repulsion ^ in the 
former, while the latter corresponds to completely uncorrelated energy levels. 

The system of free electrons in a random potential is a unique example where both uni- 
versal statistics can exist depending on the strength of disorder. For low values of impurity 
concentration rii < n c or energies E above the mobility edge E c , the electron eigenstates are 
extended and overlap well with each other. The spectral correlations in such a metallic phase 
is well described by WD statistics 2,3 * ) . With disorder increasing the system of space dimen- 
sionality d > 2 exhibits the Anderson metal-insulator transition 4 ) at rij = n c or E = E c . 
Below this transition, in the insulator phase, the typical wave function is localized within the 
volume £ d which can be centered in an arbitrary point r# throughout the sample of size L ^> £. 



Since two neighbouring in energy eigenstates are typically localized in different points in space 
\ r E — r E' \ ~ L, there is almost no overlapping between them. As a result, the energy levels are 
almost uncorrelated and obey the Poisson statistics. 

There is, however, a critical region \E — E c \ < SE near the mobility edge E c where £ oc 
\E — E c \~ u > L. The number of eigenstates in this region is equal to N c = (SE)/A oc L d ~ x l v ', 
where A = (pL^)" 1 is the mean level spacing, and p is the mean density of states. According 
to the well known Harris criterion vd > 2. Therefore the number of levels N c can be arbitrary 
large in the thermodynamic (TD) limit L — > oo, and the notion of critical level statistics is 
thus meaningful. 

A physical argument that the statistics of eigenvalues is determined by the statistical prop- 
erties of overlapping eigenstates, implies that the critical spectral statistics at the Anderson 
transition might be a representative of the new broad class of universal spectral statistics which 
is different both from WD and from Poisson statistics. The origin for this statement is in the 
special space structure of critical eigenstates which posess the multifractality 5 ~ 7 \ 

Very roughly, one can imagine such a multifractal structure by considering the space regions 
Q where the modulus of a typical wave function |\l/(r)| is bigger than a given number M. Then 
for not too large M the space pattern of Q will be qualitatively different in the critical spectral 
region as compared to metal and insulator. 

In metal, Q covers the whole sample volume L d which is "filled" densely by the extended 
wave function. In the insulator, Q is only a small fraction (t;/L) d of the sample volume that 
tends to zero in the TD limit. In contrast to both these typical pattern of f2 for a 

critical eigenstate is a sparse fractal cluster which is spread throughout the sample but whose 
total volume is proportional to L d * <C L d . "Multifractality" corresponds to the dependence of 
the cluster fractal dimension d* < d on the " cut-off level" M. 

More rigorous definition of the multifractality can be done in terms of the L-dependence of 
the moments of inverse participation ratio: 



E p (E)=^2(\^n(r)\ 2p 5(E-e n ))^< 



£-<2(p-i). metal- 
const.; insulator; 
k L-^(p)(p-i) ; critical. 



(1) 



where e n is an exact eigenvalue, (.) denotes the ensemble average, and p is an integer. 

As at large p big slashes of the wave function are enhanced, the dependence of the fractal 
dimension d*(p) has the same nature as the dependence d*(M) in the above qualitative picture. 

The universality of the critical level statistics first mentioned in 8 ) arises from the scaling 
picture of the Anderson transition. The necessary element of any such picture is an existence 
of a universal fixed point in the space of relevant parameters which determines all physics in 
the critical region. Within the one parameter scaling by Abrahams, Anderson, Licciardello and 
Ramakrishnan 9 \ the only relevant parameter is the dimensionless conductance g(L). If the 
short-scale dimensionless conductance g(l) (where I is the elastic scattering mean free path) 
is smaller than the conductance at the fixed point g c , the system flows (with increasing L) to 
the insulating state, g(L) — > 0, and thus the spectral statistics in the TD limit approaches 
Poisson statistics. If g(l) > g c , the system flows to the metal state, g(L) — > oo, where the 
spectral statistics is WD. However, as long as we are in the critical region, the dimensionless 
conductance is unchanged g(L) m g c , and thus the critical statistics is universal and does not 
depend on the system size L. 

Note that the spectral statistics take a universal form only in the TD limit when such 
parameters as l/L vanish. Since the mean level spacing A vanishes in this limit too, the 
meaningful description of the level statistics can be done only if all energies are measured in 
units of A, so that separate levels can be resolved, no matter how small is A. 



Physics of critical phenomena which exhibit scaling is described by a number of power-law 
dependences characterizing by the corresponding critical exponents. This is true for the scaling 
theory of the Anderson transition too. Within the simple one-parameter scaling 9 \ there is only 
one non-trivial critical exponent v which determines the localization length £ oc \E — E c \~ u . 
All other critical exponents can be expressed in terms of z/, the symmetry parameter (5 and the 
dimensionality of space d. 

In the recent works 10,11 we studied analytically the behavior of the correlation function 
R{s) of two levels separated by an energy interval uo = sA using the simple one-parameter 
scaling. For this correlation function it is the power-law tail R{s) = — As~^ at s ^> 1 which 
proves to depend on properties of the fixed point g = g c . The new "spectral" critical exponent 
p has been related in 10,11 ) with the critical exponent v as follows p = 1 + {yd) -1 . 

Multifractality of wave functions is described by a whole set of additional critical exponents 
r) p — d — d*(p) which were set zero in the works 10 > n ). The main objective of the present 
contribution is to show where and how they can show up themselves in the spectral statistics. 



II. LEVEL NUMBER VARIANCE AND THE NORMALIZATION SUM RULE. 

The first statement 12 -* on the critical level statistics made in 1988 was about the behavior 
of the level number variance £ 2 (A0 = {(N — N) 2 ) in an energy window of width NA (1 <C 
N <C N c ) centered at the mobility edge E = E c . This quantity is an integral characteristic of 
two-level correlations that can be related 1 ' 10 ^ to R(s) as follows: 

£ 2 (iV) = ((N - AO 2 ) = J__ (N - | a |) R(s) ds, (2) 

where the two-level correlation function R(s) is defined in terms of the exact density of states 
p{E): 

R{s)p 2 = (p(E + sA)p(E)) - (p(E + sA))(p(E)) (3) 

Making use of the perturbative expression for R(s) derived in 13 ^ which is valid in metal for 
g(L) ^> 1, and then assuming its validity up to the fixed point g(L) — g c ~ 1, it was argued 
in 12 ) that the level number variance ^(N) = aN should be linear in with the coefficient 
a < 1. Since then the linear dependence of £2 (AO for A" ^> 1 has been confirmed by a number 
of numerical simulations 12 > 14 > 15 on the tight-binding Anderson model. 

Thus the S 2 statistic at the mobility edge seems to be similar to that in the insulator where 
according to the Poisson law S 2 (A^) = N, and differs drastically from the S 2 statistic in metal 
where RMT predicts ^ S 2 (A0 oc lnA^. It is in contrast to other statistics such as the two-level 
correlation function 10 > n > 16 ) or spacing distribution function 8 ' 14 ' 15 ) which in the critical region 
behave like s 13 at small s exhibiting the level repulsion similar to that in metal. 

A challenging issue about the critical level statistics is that two levels repel each other at 
all scales, yet a big number of levels in an energy window fluctuates as if a fraction of them, 
aN, are completely uncorrelated. 

An attempt to understand deeper the nature of the linear term in S 2 statistic leads us to 
re-examening of the normalization sum rule. Suppose we have a finite sample with the total 
number of degrees of freedom (and the total number of eigenstates) equals / p(E) dE = Af . 
In contrast to the density of states p(E) that takes different values for different realizations of 



disorder, the number N does not fluctuate at all. Using this fact and integrating Eq.(||) over s 
we immediately obtain: 



+00 

R M {s) ds oc (Afp(E)} - (AO (p(E)) = 0, (4) 

-00 

where the subscript M implies that we consider a sample of a finite size L. 
On the other hand, differentiating Eq.(^) with respect to N we have: 

Looking at Eqs.(Q) and one could draw a conclusion 10 ) that dT, 2 /dN — > as N — > 00 and 
thus the linear term in £2 is excluded. This statement is in fact wrong. The point is that 
studying the critical level statistics we are bound to fulfil the inequality N < N c <^ Af. Thus 
at fixed L we cannot take the limit N — > 00. There is always the remainder integral from 
|s| = A^ c to |s| ~ A/" in Eq.([|) where the integrand, Rj^(s), is neither critical nor universal. 
This difficulty could seem to be circumvented if we take the TD limit as we described before. 
Then both N c and M go to infinity, and the remainder integral should vanish provided that 
the function Rj^f(s) decreases rapidly enough with increasing s. It is the last condition that 
makes the situation non-trivial. If the function R_\f(s) has a small but slowly decreasing tail 
R_m(s) ~ Af k ~ 1 s~ k (0 < k < 1), it becomes possible 17 ' 18 ) that in the TD limit R_\f(s) vanishes 
on a segment \s\ e (iV c — > 00, M — > 00) but the integral of it is not. Moreover, since the function 
Rn(s) becomes universal in the TD limit in which N c /Af — > 0, the integral Rat{s) ds is 
universal too and must be exactly compensated by the remainder integral in Eq.(^). This 
means that the TD limit of the remainder integral is universal 17>18 ) despite the non-universal 
integrand. 

One could avoid the above long description, just saying that taking the TD limit and doing 
the integral in Eq.([|) do not necessarily commute. Having an objective to study the universal 
spectral statistics one has to take the TD limit which eliminates the non-universal long tails 
in Rj^-(s). Then it is a non-trivial question if the normalization sum rule proved rigorously for 
finite J\f survive taking the TD limit. It does not if the remainder integral tends to a non-zero 
limit —a as L —>■ 00. This is just the deficiency of the sum rule: 

+00 

Roo(s) ds = a. (6) 

-00 

Note that in doing the integral in Eq.(^) there is no problem with taking the TD limit, since 
the limits of integration are finite and fixed. Now, using Eq.(^) and assuming the power-law 
tail in Roo(s) = —As^^ with 2 > fi > 1 we arrive at: 

^ { N) = aN + ¥ - M - W) +mnSt . (7) 

The second term in Eq.(^) is totally determined by the form of the tail in the two-level 
correlation function and reduces to IniV in the RMT limit \i — > 2. However, the linear in N 
term is due to the violation of the normalization sum rule in the TD limit. 



III. THE SPECTRAL FORM-FACTOR AND THE RETURN PROBABILITY. 



A convenient equivalent formulation of the sum rule violation Eq.@ can be done in terms 
of the t — > limit of the spectral form-factor 2nK(t) = j_™ Roo(s) exp[— ist] ds: 

a = 2ttK(0). (8) 

The latter can be related to the probability L~ d p(t) for a particle to return to the initial point 
for the time T = th/A. A perturbative relationship 19 ) between Kit) and p(t) which is valid 
in metals for t < 1 ( or T <C h/A) and is equivalent to the diagrammatic approach 13 \ reads: 

K(t) = (2ny 2 \t\p(\t\), (9) 

where p(t) is given by the diffusion propagator -P(q, oS): 

1 f 2 c dt i 

Pit) = j^- 2 J ^(Gf +sA (r,r)G £ A (r,r))e-^ = -j — £P(q,u,) exp[-iut/A]. (10) 

Here G^ A ^ (r, r') are exact retarded (advanced) Green functions of electrons in a random po- 
tential, and the factor 2/(3 accounts for the number of 'massless' diffusion modes for a given 
symmetry class. 

If the multifractality of wave functions is not taken into account, then Eq.(^|) is sufficient to 
describe quantitatively the long-range part of two-level correlations in metal and to establish 
the power law tail in R(s) up to the numerical pre-factor in the critical region. Indeed, in metal 
we have P(q, us) = {Dq 2 — ioo)' 1 , where the momentum q is quantized qi = jni (rii = 0, 1, 2...). 
For sufficiently small u D/L 2 , or a time which is much larger than the diffusion time, the 
diffusing wave packet initially <5(r)-shaped at the origin, is spead homogeniously throughout 
the sample. It is equivalent to considering only the zero-mode with q = in the sum, Eq. (|10|) . 
Then one immediately arrives at p(t) = 2/(3 and K(t) = \t\/2n 2 (3. This corresponds to the 
WD tail in the two-level correlation function R(s) = —l/ir 2 j3s 2 . For smaller times T L 2 /D, 
the wave packet is spreading diffusively, and p(t) oc {Dt)~ d l 2 . The diffusive propagation of a 
wave packet corresponds 19 - 1 to the Altshuler-Shklovskii regime of level correlations 13 - 1 where 
R(s) ocs-( 2 ~ d / 2 l 

It turns out 10 ' 11 - ) that near the Anderson transition, where the correlation/localization 
length is large £ 3> /, one more regime is possible, which is absent both in a good metal and in 
a strongly localized state. It corresponds to diffusion at a length scale r £. According to the 
one-parameter scaling 4 ) , diffusion in this regime is anomalous with the scale-dependent diffusion 
coefficient D{r) = p~ l r 2 ~ d g{r) related to the the scale dependent conductance g{r) = g c + 5g(r), 
where: 

Sg(r) ~ 9c(r/0 1/u - (11) 

If one sets g(r) = g c as the first approximation, then one arrives at the anomalous diffusion 
where r d oc t is proportional to time for T < h/A^ = hp^ d . Then the quantity tp{t) oc t/r d 
reduces to a constant which is irrelevant for the power-law tail in R(s). For s 3> A^/A (which 
corresponds to t A/A^ ) this tail turns out 10 ' n ) to be determined by the correction Sg(r), 
Eq.(0): 



K(t)~tp(t) = -L ~ I + const. R(s) ~ I .-0+A>, (12) 

g(r) g c \ A J g c \ A J 



where 



A \l; 



(13) 



IV. THE MULTIFRACTALITY OF WAVE FUNCTIONS AND SCALING. 



However, the relationship Eq.@ does not give even a qualitative description of the critical 
level statistics if multifractality is taken into account. In order to clarify this point let us consider 
the representation of R(s) in terms of diagrams with exact diffusion propagators P(q, ui) and 
exact vertex parts r 2n +i({qfc}; oS) which depend on 2n — 1 momenta qfc. A typical {2n — l)-loop 
diagram is shown in Fig.l. 



P(q,,<o) 




P(-q -q -...-q ,co) 

1 1 1 2 1 zn-1 

FIG. 1. A skeleton diagram for R(s). 



The relationship between K(t) and p(t) given by Eq.(^) exactly corresponds to the one-loop 
diagram first considered in 13 ). 

It is clear that in the critical region all diagrams should be considered. Why then we believe 
that the power-law tail is correctly described (up to a numerical pre-factor) by Eq.(|9])? The 
reason is the scaling behavior of the diffusion propagator. Using the scaling of the diffusion 
constant D(r) oc r 2 ~ d one concludes that Dq 2 scales as r~ d or q . The scaling dimension of a 
vertex is equal to zero. Each loop in diagrams of Fig.l adds an extra diffusion propagator that 
scales like q~ d . This scaling factor is exactly compensated by a factor q d which arises because 
of an extra integration over momenta. Thus all diagrams shown in Fig.l have the same scaling 
dimension and should be proportional to the same power of the quantity L w = (—iptu)~ 1 / d (its 
absolute value is equal to the mean root square of the diffusive displacement for the time 1/uS) 
provided that it is the only relevant length in the system in the critical region. In case of no 
multifractality the momentum integration in all diagrams is convergent at large momenta and 
the above assumption is correct. However the situation changes if multifractality is taken into 
account. We will see below that in this case all diagrams are divergent, and the inverse upper 
momentum cut-off A -1 ~ / becomes a relevant length too. 

Let us introduce the generalized diffusion coefficient D(q,u) which obeys the above scaling: 

D(q,u;) = g c p- 1 q d ' 2 F(qL UJ ), (14) 

where F(qL u ) is a universal scaling function. A usual assumption is that this function is such 
that D(q,u>) is independent of q at qL u <C 1 and it is independent of u in the opposite limit. 
It immediately leads to D(q,ou) oc q d ~ 2 at large q. 

However, the multifractality of wave functions should exhibit itself in the large-g limit of 
the diffusion propagator P(q, u) = [q 2 D(q,uj) — ioS\~ l . Thus it is natural to assume 20 ) that in 
case of multifractality D(q,ui) oc q d *~ 2 at qL w ^> 1. Then Eq . fll"4|) leads to: 

q 2 D{q,u) = Cg cP - l q d ^L- u \ (15) 

where C is a numerical factor, and 77 = d — d*(2). For a diffusion propagator that is expressed 
in terms of the product l 1 !/™! 2 ^™! 2 (see Eq.fllOD), the relevant fractal dimension is d* = d*(2). 

An immediate consequence of Eq.(|i"5|) is that the return probability given by Eq. fllOD is 
divergent at large momenta: 



pit) ~ • ( 16 ) 

This divergency, in turn, results in the divergency of K(t) in Eq.(||). It is instructive to derive 
how this divergency arises directly from the one- loop diagram in Fig.l. In order to do that we 
note that there is a Ward identity that relates the vertex part r 3 (q, o>) with the derivative of 
the diffusion propagator with respect to lo: 

r 8 (q, W ) = P- 2 (q,-) = --^[P-\ q ,u)]. (17) 

Taking into account Eq.flTTD and the ^-dependence of the generalized diffusion coefficient at 
qL^ ^> 1 which follows from Eq . (|15|) , we arrive at: 



r 3 (q,u;)=r^ +iC-^g c {qL») 1 



where I3 = i is a constant. 



Eq.(P"8"D demonstrates the generic property of vertex parts r 2n+1 ({q}, u>) at qL u 1: for 
r) all of them are growing functions of q. It is the second term in Eq.flTjp that leads to the 
divergent one-loop diagram for R(s) at r\ ^ 0. 

One can continue the power-counting for an arbitrary (2n — l)-loop diagram using the 
recursive relation: 



Ar 2 n + i({q},q'q»P(q» = A nP T 2n ^({q},u;), (19) 
where A n is a numerical coefficient. Thus we arrive at: 

[p(s)}-\ (20) 



T2n+1 



1 2n-l 

where 



p{s) = p- 1 J d d qP{q,uj) = "J dtp{t)e ist ~ {AL) r 's~^. (21) 

Now one can estimate an arbitrary (2n — l)-loop diagram for R(s). All of them are divergent 
at large momenta so that we consider only the second term in Eq.(^) to solve the recursive 
relation Eq. (|20|) . Using also Eqs.(|l5|) and (^l|) we have: 



#2n-i(s) 



ft / d"q*P(q*) 

k=i 



1 



5(qi + ... + q 2 n)r 2 „ + i({q}; u) ~ - [p(s)} n . (22) 



s 



The main lesson 21 ) we learn from the above crude power-counting is that as soon as i] ^ 0, 
the power of the divergent parameter Q = (AL) V increases with the number of loops 2n — 1 in 
a diagram of Fig.l. Were it not so, there would be no way to get a finite TD limit of R(s) as 
(AL) — > oo. However, Eq.(E2|) tells us that the possible way out is to assume 21 ) that the sum 
R( s ) — 2^Li^2n-i( s ) is ver y much like the geometric progression Q + Q 2 + ... = Q/(l — Q) 
which is finite in the limit Q — > oo. 

Another lesson is that within the power counting the two-level correlation function R{s) 
can be represented as a series in p(s) which is the Fourier-transform of the return probability 
pit), all the divergency being absorbed in p(s) or p(t). Thus the possible expression of K{t) 
in terms of p(t) that realizes the idea of the 'geometric progression' should be of the form 
K(t) oc \t\p(t)/(l + L\p(t)]), where L\p(t)] is a linear in p(t) operator. 



V. AN IMPROVED RELATIONSHIP BETWEEN K(T) AND P(T) AND THE SUM 

RULE VIOLATION. 



An expression of this kind has been recently derived 22 - ) using the idea of parametric diffusion 
in the eigenvalue space caused by a random perturbation. It reads: 

m = (2Kylll\lp(t>)dt>- (23) 

By the derivation, Eq.(^) applies only for small times T <C h/A (t <C 1). However, in contrast 
to Eq.(^), it gives a qualitatively correct constant limit lim t ^ 00 i ; r(t) = 1/4-7T, which is yet a 
factor of two smaller than the exact value K ex (oo) = 1/2tt. 

Now we make use of Eqs.(^) and flTB]) to relate the deficiency of the sum rule, Eq.(|8|) with 
the multifractality exponent rj. Implying the TD limit and the fact that p(t) is divergent in 
this limit, we retain in the denominator of Eq.(^) only the term proportional to a divergent 
parameter Q = (AL) ri . Then we immediately arrive at 26 ^: 

(24, 

This expression is a central in our consideration. It gives an answer to a question why the 
normalization sum rule is violated in the critical region. According to Eq. (|24]) it is only because 
of the multifractality of wave functions near the mobility edge. 

From the definition, Eq.(|l|), of the fractal dimensions d*(p) one can conclude that the WD 
statistics in metals correspond formally to d*(p) = d and the Poisson statistics in insulator 
correspond to d*(p) = 0. Indeed, in case of WD statistics, the normalization sum rule, Eq.(|), 
is known to survive the TD limit and no linear term arises in the level number variance, Eq.(^), 
as it follows from Eq.fl2|) at d *( 2 ) = d - 

However, Eq.([24"t) is not able to reproduce a correct result a = 1 in case of Poisson statistics. 
The reason is that there was an assumption made in the derivation of Eq.(p3|) that implies a 
sort of decoupling between statistics of eigenvalues and eigenfunctions. Such an assumption 
can be correct only for relatively weak space structure of eigenfunctions, that is for small values 
of rj and away from the strong localization regime. Thus Eq. ( p4|) can be considered as the first 
term in the expansion of a in powers of the small parameter r)/d. 

The validity of Eq.fl24|) can be checked using the exact results obtained in 2D systems in 
the regime of weak localization L -C £ = le 90 , where go = (2n 2 pD) 3> 1. It is known 5 - ) that 
in this regime the moments of inverse participation ratio behave in the same way as for the 
critical eigenstates in Eq. ([!]). The corresponding spectrum of fractal dimensions has been found 
recently in the framework of the nonlinear supersymmetric sigma model 23 ^ and coincides with 
the earlier renormalization group result 5 ' 24 ^: 

cT(p)=2--£-. (25) 

On the other hand, as long as the parameter g^ 1 is small, there is no problem to calculate the 
behavior of the two-level correlation function Rjy(s) in the region of large energies u = sA 3> 
D/L 2 (which are however smaller than the Fermi energy ef 3> h/r). It has the form 25 ' 17 ); 



471^0 



s 2 5 2 

(26) 



l + s 2 5 2 

where the parameter 5 = (rA/h) is proportional to A/" -1 and vanishes in the TD limit 



The expression, Eq.fl2"6|) is valid in the ballistic region (a; ~ ^/ r ) too. It illustrates the 
existence of a small but slowly decreasing tail which has been discussed in Section II. For 
s < ~ A/", the tail is proportional to 5ln(5~ 1 /s) and is almost s-independent, while for 
s ^> 1/5 it decreases rapidly. This tail vanishes in the TD limit but the integral of it tends to 
a finite limit for any fixed lower cut-off a: 



r°° I 
lim / R M ( s )ds = -a = - — . (27) 



Thus for the 2D system we know (to the leading order in g x ) expressions both for d*{p) and 
a calculated independently. It is easy to check that the relationship between them is really the 



one that follows from Eq.(24). 



VI. DISCUSSION AND OPEN QUESTIONS. 

In the above consideration we have shown that the quasi-Poisson linear term in the critical 
level number variance ^{N) arises because of the multifractal space structure of wave functions 
near the mobility edge. Thus we conclude that one of the simplest form all eigenvalue statistics, 
the level number variance, contains an information about the fractal dimension d*{2) = d — rj 
which characterizes the non-trivial statistics of eigenf unctions. 



It is interesting to compare the value of r\ predicted by Eq. (|24[) with that found by the direct 
numerical investigation 27 ) of statistics of critical eigenfunctions in the 3D Anderson model. In 
order to do that we use the numerically obtained slope a in the E 2 (iV) statistic at the mobility 
edge. Simulations on the 3D Anderson model 12 > 14 > 15 ) give the value of a ~ 0.25 — 0.30 with 
the best estimation 15 - ) being a = 0.27. Then Eq.([24j) predicts r] = 1.62 ± 0.15. This is in 
a surprisingly good agreement with an old result d*(2) = 1.7 ± 0.3 of the direct numerical 
evaluation 2? ) of d*(2). 

So far we were interested only in a limit a = 27iTimt^o K(t) in Eq. (|23"D . Within the scaling 
approach we used in Section V, we were able to find only leading power-law dependences on 
t. For that reason, the spectral form factor K(t) was independent of t (in the TD limit) 
whatsoever. Small t-dependent term in 2isK.it) = a + would lead to a power-law tail in 
R{s) oc s~ M in the same way as in Eq.(|l2|). A natural question which arises in this connection 
is whether the result, Eq.(|T2|) survives if multifractality and the improved relationship Eq.(|2BJ) 
between K{t) and p(t) is taken into account instead of Kit) oc \t\p(t). 

In order to answer this question let us substitute g c for g{r) = g c + 5g(r) in Eq. fllBD , where 
5g(r) is given by Eq.(|ll|) with r d cx t. The result of such a substitution is the multiplication of 
pit) in Eq.(||) by the same factor of [1 + A(A^/A)^] as in Eq . (|12|) . Now we use Eq. (|23|) to 
arrive at the result for K{i) that differs only in the numerical coefficient in front of t^ from 
Eq. (|T2|) . Thus the multifractality and the improved relationship between K(t) and p(t) do leave 
the result Eq . (|12"D unchanged. 

However, there could be additional power-law corrections to Eq. (|T6|) other than that coming 
from the correction to conductance 5g(r). A possible source of such corrections is the scaling 
function in Eq.(|TJ|). For L = £ = oo it depends only on qL^ but at finite L and £ it has 1/L or 
l/£ corrections. In particular for qL^ ^> 1 it should cross over from (qL u )~ v for L w L < £ 
to (qL)'* 1 for £ ^> L w ^> L. One can imagine a cross-over of the type F(gL w , L u / L) = 
(qL u )-^[l + (L U /L)<Y>K. It is characterized by a new cross-over exponent ( which play a role 
similar to \ jv in corrections to conductance considered above. Thus the question about the 
exponent /i in the power-law tail of R{s) is not that simple as it seemed earlier and requires 
further investigations. 
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Statistique spectrale a la transition d'Anderson: Multifractalite des fonctions 
d'ondes et violation d'une regie de somme de normalisation. 



La statistique des niveaux d'energie d'electrons dans un potentiel aleatoire est consideree dans une 
fenetre d'energie critique autours du bord de mobilite. On montre que la multifractalite des fonctions 
d'onde critiques a pour consequence la violation d'une regie de somme de normalisation dans la limite 
thermodynamique et donne un contribution quasi-poissonien < (5(N) 2 >= a < N > a la variance du 
nombre de niveaux. La deficience a = rj/2d a la regie de somme est reliee a l'exposant multifractal 
rj = d — d* (2) dans un espace d-dimensionnel. 



